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THEORIES AND HEAT PULSE EXPERIMENTS OF 
NON-FOURIER HEAT CONDUCTION 

P. VAN 


Abstract. The experimental basis and theoretical background of non-Fourier 
heat conduction is shortly reviewed from the point of view of non-equilibrium 
thermodynamics. The performance of different theories is compared in case of 
heat pulse experiments. 


1. Introduction 

There are several extensive reviews of heat conduction beyond Fourier [Dill 
111 1H. When non-equilibrium thermodynamics is concerned, usually Extended 
Thermodynamics is in the focus of surveys in this subject. In this short paper heat 
conduction phenomena is reviewed from a broader perspective of non-equilibrium 
thermodynamics. The performance of different theories is compared on the example 
of heat pulse experiments. 

There are two preliminary remarks to underline the particular point of view of 
this work. The first remark concerns the so called heat conduction paradox of the 
infinite speed of signal propagation in Fourier theory. Infinite propagation speeds 
are frequently mentioned as a disadvantage mu and this is a main motivation for 
looking symmetric hyperbolic evolution equations in general |9]. However 

- the relativistic formulation of heat conduction is a different matter and 
parabolic equations may be compatible with relativity [Hiniin]. 

- The finite characteristic speeds of hyperbolic equations are determined by 
material parameters. Therefore, in principle, in nonrelativistic spacetime 
they can be larger than the speed of light. 

- In case of parabolic equations the speed of real signal propagation is finite, 
e.g. because the validity condition of a continuum theory determines a 
propagation speed. Moreover, the sensitivity of the experimental devices is 
finite, therefore they can measure only finite speeds muni mug. 

Our second remark is about the role of non-equilibrium thermodynamics as a 
frame theory of several physical disciplines. Here the fundamental assumption is 
the validity of the second law of thermodynamics in the form of an entropy balance 
with nonnegative production. If microscopic or mesoscopic theories do not want 
to violate macroscopic second law, they must suit the requirements derived from 
the phenomenological conditions. This is the extension of the classical universality 
of the absolute temperature. The validity of the second law is a weak assumption 
from a physical point of view but introduces remarkably strict restrictions for the 
evolution equations. 
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According to these arguments it is no unreasonable to look for parabolic or at 
least non-hyperbolic extensions of the Fourier equation of heat conduction that 
fulfill basic principles of continuum physics. 

In the following section we survey the most important experimental observations 
of non-Fourier heat conduction. Then the relevant kinetic and thermodynamic 
theories are surveyed. Their assumptions, properties and modeling capabilities are 
compared on the example of heat pulse experiments. 

2. Experiments 

2.1. Helium II. The hrst experimental observation of non-Fourier heat conduction 
was the measurement of wave like propagation of heat in liquid Helium II by Peshkov 
m- The experiment was motivated by two-fluid theories of Tisza m and Landau 
[18]. Landau suggested that wave like propagation is a property of phonon gas 
and introduced the terminology ’’second sound”. The works of Cattaneo, Morse- 
Feshbach and Vernotte [3111110] indicated that the related phenomena may be 
more general, there is an ’’inertia of heat”. 

2.2. Low temperature solids. Phonon gas theory predicted the existence of sec¬ 
ond sound in solids. The key aspects of the observation were the modeling of 
phonon scattering by kinetic theory midi HO] and the application of heat pulses. 
The derivation of a dissipative extension of the Maxwell-Cattaneo-Vernotte equa¬ 
tion from kinetic theory, the Guyer-Krumhansl equation [22] . provided the ’’window 
condition”, when dissipation is minimal. A careful preparation of crystals with par¬ 
ticular properties resulted in observation of second sound in solid ^He, "^He, NaF 
an Bi crystals m m m dz]. In some experiments second sound and ballistic 
propagation - heat pulses propagating with the speed of sound - were observed 
together [28]. At lower temperatures ballistic propagation appears without second 
sound. 

The correct modeling of the observed parallel ballistic, wave like and diffusive 
propagation of heat in a uniform theoretical framework is the most important 
benchmark of the theories. 

2.3. Heterogeneous materials at room temperature. The low temperature 
heat pulse measurements in dielectric crystals exploit well understood microscopic 
mechanisms. At room temperature various dissipative processes suppress these 
effects. However, some experiments with heterogeneous materials indicated the 
possibility of non-Fourier heat conduction at room temperature, too [13[30]. These 
experiments were not confirmed, more properly the attempts of exact reproduction 
of these experiments are contradictory [SlEllESlIll]. 

2.4. Small size. Nano heat conduction is a popular research field with interesting 
new observations [35] . Non-Fourier effects are in principle enhanced by reducing the 
size of the samples and the speed of the phenomena (36] |37| [38] . Guyer-Krumhansl 
equation is extensively analyzed in this respect (39] [40] . 

2.5. Relativistic experiments. Relativistic fluid effects became experimentally 
available in quark-gluon plasma. Experiments in RHIC and LHG confirmed the 
existence of a dissipative relativistic fluid [41]. In relativistic theories the inertial 
effects cannot be neglected, viable theories incorporate wave like heat propagation, 
due to stability problems gaiisiiii]. However, in case of ultrarelativistic speeds 
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the flow is energy dominated, therefore it is reasonable to use Landau-Lifshitz flow- 
frame, where propagation of energy defines the flow. This choice eliminates the 
possibility of heat conduction and only viscous dissipation is possible. According 
to thermodynamic arguments, the choice of flow-frames is not arbitrary and in case 
of lower speeds heat conduction may be separated from viscous dissipation |45) . 

3 . Theories 

The three main approaches are kinetic theory, non-equilibrium thermodynamics 
and continuum mechanics combined with non-Fourier heat conduction. Up to now 
the relation of these general approaches is not clear. Fortunately these general 
frame theories offer several possibilities of comparison, validation and extension far 
beyond heat conduction phenomena. In particular we will see, that the performance 
of these theories in modeling ballistic and diffusive heat propagation is different. 

4. Kinetic theory 

As it was already mentioned, kinetic theory played an important predictive role 
in understanding and designing proper experiments with the help of identified mi¬ 
croscopic heat propagation mechanisms. In particular the relation of pure kinetic 
and continuum approaches is important in this respect. Kinetic theory introduces 
a hierarchical structure of macroscopic field quantities with coupled balance form 
evolution equations. In this balance structure current density at a given level is a 
source density at the next level. The tensorial order of the n-th variable is n. For 
phonons with the Callaway collision integral one may obtain the following system 
of equations in one spatial dimension ([9] p349): 

n = 0, 
n = 1, 

n >2. 

Here dt and dx are the time and space derivatives, Un is related to the nth 
momentum of the one particle probability distribution function by constant multi¬ 
pliers, T/j and Tjv are relaxation times, c is the Debye speed. A truncated hierarchy 
at the third moment, = 0 of n > 3, results in the following set of equations: 

(2) dtUQ + cdxUi = 

( 3 ) dtui + ^cdxUQ + cdxU2 = 

4 

(4) dtU2 + —cdxUi = 

15 

The macroscopic fields are introduced by the following definitions: the energy 
density e = Hcuq, heat flux q = hc^ui and the moment of the heat flux Q = hc^ui. 
One may apply the (approximate) caloric equation of state e = pcT with con¬ 
stant density p and specific heat c. Moreover, it is convenient to redefine the 
coefficients by introducing the g-relaxation time Tq = r^, the Q-relaxation time 
tq = XR+TN ’ Fourier heat conduction coefficient A = pctrc^ and two ad¬ 
ditional phenomenological coefficients k2i = trc^ and ki2 = 4 t / 15 . Then the 
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previous equations may be written as: 


(5) 

dte dxQ 

(6) 

Tgdtq -b XdxT -b K 2 AQ 

(7) 

TgdtQ -b Ki2dxq 


If TQ = 0, one may obtain the Guyer-Krumhansl equation eliminating Q from 
([6])-([7]). If K 2 i = 0 then ([6]) is the Maxwell-Cattaneo-Vernotte equation. 

In Rational Extended Thermodynamics (RET) the coefficients of the above set 
of held equations are calculated exactly as it was shown above. Remarkable that 
one needs only three parameters, the Debye speed and the R and N process related 
relaxation times [9l |46l |6] . 

Extended Irreversible Thermodynamics (EIT) is a related but different approach, 
where the structure of the equations is important, but the values of the coefficients 
are considered mostly undetermined, therefore e.g. in ©-(O the number of pa¬ 
rameters is higher liEzl 1^ . In RET the microstructure is unavoidable in EIT 
the microstructure is hexible. In RET the coefficients are to be calculated from a 
microscopic model, in EIT some coefficients are to be measured. RET is a local the¬ 
ory by construction, EIT may be weakly nonlocal. RET is strictly compatible with 
kinetic theory, EIT is weakly compatible with kinetic theory. In EIT the derivation 
of the Guyer-Krumhansl equation either refers to the momentum hierarchy [47j or 
weakly nonlocal extensions in a pure non-equilibrium thermodynamic framework 

[nmg. 

The experimental results of heat conduction of mixed ballistic-wavy-diffusive 
propagation can be modeled in this framework. However, with the restriction of 
using only 3 parameters RET requires about 30 (!) moments in order to obtain 
correct propagation speeds both for ballistic heat propagation and second sound 
[igE]. In EIT, considering some coefficients as phenomenological parameters, in 
principle it is possible to incorporate correct propagation speeds for the second 
sound and also for the ballistic phonons with the above set of equations. However, 
in this case one cannot stop at the level of the Guyer-Krumhansl equation, and, most 
importantly, more than 3 parameters are required. A hierarchical phenomenological 
background theory is necessary to justify these assumptions. 

5. Non-equilibrium thermodynamics 

Irreversible thermodynamics as held theory was established by Eckart, with the 
assumption of local equilibrium [501 EH Ea [53]. Maxwell-Gattaneo-Vernotte type 
extensions require a deviation from local equilibrium. In this paper non-equilibrium 
thermodynamics is a collective nomination of thermodynamic theories based on the 
direct exploitation of the entropy inequality with or without the requirement of local 
equilibrium. Classical Irreversible Thermodynamics, Extended Thermodynamics, 
Thermodynamics with Internal Variables are special theories of non-equilibrium 
thermodynamics [46] . 

The hrst basic assumption extending non-equilibrium thermodynamics beyond 
local equilibrium is that the thermodynamic fluxes of Classical Irreversible Ther¬ 
modynamics are state variables. Then Maxwell-Cattaneo-Vernotte equation is a 
consequence of the second law by thermodynamical flux-force linearization. This 
idea was first proposed by Muller [54]. One of Muller’s argument was the com¬ 
patibility with kinetic theory. An other independent way is based on the seminal 
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treatment of discrete systems by Machlup and Onsager, who introduced kinetic 
energy of thermodynamic state variables [55) . The continuum generalization of 
Gyarmati is based on variational considerations and therefore introduced thermo¬ 
dynamic fluxes as independent variables instead of time derivatives [56j . Later on 
this idea was further generalized by introducing arbitrary internal variables for the 
same purpose [ST) [58] . 

As it was mentioned previously, thermodynamic fluxes as state variables may lead 
to the Maxwell-Cattaneo-Vernotte equation, but to obtain the Guyer-Krumhansl 
equation requires further assumptions. The modification of the entropy current is 
a straightforward idea in this respect. The idea, that the entropy current density is 
not always the heat flux divided by the temperature is natural in mixtures and also 
in Extended Thermodynamics [42] . but may be natural in other generalizations, too 
[siinniini]- To obtain the non-equilibrium thermodynamic counterpart of ®-© 
it is reasonable to assume that the additional fields are the heat flux vector q and 
a second order tensorial internal variable Q and use the Nyiri form generalization 
of the entropy current density [62l [63] . 

In this paper we restrict ourselves to rigid heat conductors in a one dimensional 
treatment, therefore the time derivatives are partial and the space derivatives are 
one directional. For a more complete treatment see [64) . Our starting point is the 
balance of entropy, which results in nonnegative production when constrained by 
the balance of internal energy ®: 

( 8 ) dts + dxJ > 0 . 

Here s is the entropy density and J is the entropy current density. 

Deviation from local equilibrium will be characterized by two basic constitutive 
hypotheses: 

-- We assume that the entropy density depends quadratically on the additional 
helds: 

(9) s(e, g, Q) = Seq{e) - 

where mi and m 2 are constant, nonnegative material coefficients. The 
derivative of the local equilibrium part of the entropy function Seq by the 
internal energy is the reciprocal temperature: 

dSeq 1 

~d^ ^ T' 

The quadratic form is introduced for the sake of simplicity and may be 
considered as a first approximation. The coefficients mi and m 2 are non¬ 
negative because of the concavity of the entropy function, that is, thermo¬ 
dynamic stability. 

- We assume that the entropy flux is zero if g = 0 and Q = 0. Therefore it 
can be written in the following form: 

(10) J = bq + BQ. 

Here b is originated from a second order tensorial constitutive function and 
B from a third order one. They are the current multipliers [62) . 
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Here the basic fields are T, q and Q, the constitutive functions are b and B. The 
entropy production can be calculated accordingly; 

dts + dxJ = -^dxq - rrnqdtq - m2QdtQ + bdxq + 
qdxb + BdxQ + QdxB 

(11) = (^b - q+{dxb - midtq)q-{dxB - m 2 dtQ)Q + Bda;Q >0. 

In the last row the first and third terms are products of second order tensors, the 
second term is a product of vectors and the last term is of third order ones. It is 
not apparent in our one dimensional simple notation. The time derivatives of the 
state variables q and Q represent their evolution equations, which are considered 
as constitutive quantities together with the current multipliers b and B. Therefore, 
one can identify four thermodynamic forces and currents in the above expression 
and assume linear relationship between them in order to obtain the solution of the 
entropy inequality. In case of isotropic materials only the second order tensors can 
show cross effects, the vectorial and third order tensorial terms are independent. 

The linear relations between the thermodynamic fluxes and forces result in the 
following constitutive equations: 


( 12 ) 

midtq - dxb = 

-hq, 

( 13 ) 

m2dtQ - d^B = 

-kiQ + ki2dxq, 

( 14 ) 

b-— = 
T 

-k2iQ - 1 - k2dj:q, 

( 15 ) 

B = 

TidxQ • 


The li^ki,k 2 ,n coefficients are nonnegative and also K = ^ 1 ^ 2 —^ 12^21 > 0, because 
of the entropy inequality (ITT|) . The constitutive equations (fT^ - dTSl) together with 
the energy balance and the caloric equation of state e = pcT form a solvable set 
of equations. Moreover, in case of constant coefficients one can easily eliminate the 
current multipliers by substituting them from dTT)) - (fT51) into (fT^ - dT51) and obtain: 

(16) midtq + hq - k 2 dlq = dx^-k 2 idxQ, 

(17) rn2dtQ + kiQ - ndlQ = ki2dxq- 

Here denotes the i-th partial derivative by x. (HSli-dlTl) is identical to 
® of the previous section introducing Xq — 777 - 1 ^/ 1 , A — (/iT ) , s: 2 i — ^ 21/^15 

tq = 11 X 2 /ki and K 12 = ^ 12 /^ 1 - If we assume that the coefficients ^2 and n are zero 
we obtain exactly ®- 0 . For the correct sign of the coefficients there one cannot 
assume reciprocity here. Guyer-Krumhansl equation is a special case when n = 0, 
m 2 = 0 and Maxwell-Cattaneo-Vernotte is obtained if fc 2 = 0 and either ki 2 = 0 
of k 2 i = 0 in addition. Several other heat conduction equations are obtained in 
this framework. Jeffreys type is included if n = 0, mi = 0, fc 2 = 0 and either 
ki 2 = 0 of k 2 i = 0. A Cahn-Hilliard type equation is derived if n = 0, mi = 0 and 
m 2 = 0 [65j . With suitable constitutive assumptions several nonlinear equations fit 
in this general framework like the thermomass transport equation [66] , the nonlinear 
Maxwell-Cattaneo-Vernotte equation derived in [39] and also the nonlinear Guyer- 
Krumhansl equation derived in m- 
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It is remarkable that the Green-Naghdi model III and II cannot be obtained here, 
however it is a valid special case if a general vectorial internal variable is introduced 
instead of the heat flux [68)1. 

6. Rational thermomechanics and some other theories 

The common property of this group is the lack of compatibility with kinetic 
theory. The rational approaches are rigorous from mathematical point of view, 
require Noll type material frame indifference and insist to classical entropy current 
density as heat flux divided by temperature. 

6.1. Second viscosity. Both ballistic propagation and second sound can be rea¬ 
sonably modeled with the help of second viscosity [Miiiniiii]- Second viscosity may 
have an imaginary part and its origin looks like essentially a kind of internal variable 
theory developed on the basis of thermostatics before the advent of non-equilibrium 
thermodynamics. The idea and the method is originally from Mandelstam and 
Leontovitch m and is well described in m- 

6.2. Jeffreys type equation. Jeffreys type equation was first suggested by Joseph 
and Preziosi [I] assuming a delayed transport in heat conduction. Later on this 
constitutive equation become very popular, because it improves the properties of 
Maxwell-Cattaneo-Vernotte equation from many point of view and because it can 
be obtained from simple (but sometimes unacceptable) assumptions [TH |351 [7S] . 

For one dimensional heat pulse experiments the Jeffreys type equation is identical 
with the Guyer-Krumhansl equation. 

6.3. Green-Naghdi equations. Green and Naghdi derived a particular constitu¬ 
tive equation for the heat flux. They have introduced an internal variable, called 
thermal displacement rate, with the particular interpretation being the time deriv¬ 
ative of the temperature m- The corresponding equations coupled to the momen¬ 
tum balance are modeling well ballistic propagation and second sound [zziiziiza 
[80]. A particular case of the model of Green and Naghdi is the existence of heat 
conduction without dissipation. 

6.4. Semi-empirical temperature of Cimmelli and Koslnski. The semi-empirical 
temperature is a scalar internal variable, in the framework of a weakly nonlocal 
theory IHI1IH1I5]- Goupled to mechanics it reproduces ballistic-wave like-diffusive 
propagation of heat [83]. However, in this theory the constitutive theory (both the 
evolution of internal variable and the heat flux - temperature relation) is postulated 
directly. 

6.5. Rational thermomechanics. There are many attempts to incorporate Maxwell- 
Gattaneo-Vernotte type heat conduction in the framework of rational thermome¬ 
chanics [mill]. It is remarkable that a nonlinear version of Maxwell-Gattaneo- 
Vernotte equation seems to fail the test of proper modeling the experiment, because 

of the lack of a certain type of dissipation, characteristic in Guyer-Krumhansl or 
Jeffreys type models [84l|49]. The theory could not meet the challenge obtaining 
compatibility with kinetic theory without the modification of the classical entropy 
current (or the kinetic theory). 
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6.6. Other approaches. There are many more methods of extending thermody¬ 
namics beyond local equilibrium. The test of these ideas usually starts with heat 
conduction. One may introduce the time derivatives of the state variables as ad¬ 
ditional state variables. However, time derivatives are frame dependent, they are 
many of them and it is not clear which one is to be used. Rigorous exploita¬ 
tion methods of the entropy principle do not help in this respect [85l [86] . Recently 
Serdyukov applied this idea for heat conduction, too [87l [88] . The ballistic-diffusive 
equation of Chen is a particular mixture of phenomenological and kinetic consider¬ 
ations [sniiso], the thermomass theory of Guo and Hou seemingly rediscovers the 
role of inertia in heat conduction |66j . 

An important class of heat conduction equations is nonlinear by construction 
[HllSaiMlISZj- Possible experimental verifications require their reliable numerical 
solution. Analytical results are crucial in this respect mm- 

These approaches are partially consistent with either kinetic theory or contin¬ 
uum mechanics. They may use or avoid the usage of rigorous entropy principle 
exploitation. Some of them, like the delayed differential equations of dual phase lag 
theory [93] failed to fulfill important thermodynamic expectations |94l|95]. None 
of them were tested with experiments of ballistic-wave like-diffusive propagation. 

7. Discussion 

In his mind provoking article Muschik writes that the reason of so many schools 
of thermodynamics is, that one may go beyond local equilibrium with different 
ways. The extension to nonlocal-nonequilibrium is not unique [96] . In this short 
survey we argued that the different extensions are not equivalent regarding their 
performance of modeling experimental results, neither regarding their theoretical 
consistency and scope. 

The challenge of non-Fouirer heat conduction is to develop a theory that is 
compatible both with kinetic theory and mechanical principles (including material 
frame indifference), passes the tests of existing experiments on parallel ballistic, 
wave like and diffusive propagation of heat and, therefore and most importantly, 
predictive in foretelling new observations and phenomena. 

This pointed review is far from being complete and the author apologizes of not 
mentioning relevant works. This is partially because of his special point of view, 
beyond his limited knowledge. 
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